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THE  THICKNESS  OF  SHOCK  AND  DETOMTION  FRONTS  AND 
ROTATIONAL  RELAXATION  AT  HIGH  TEMPERATURES 

Iqr  K,  Hansen,  D,F.  Hornig,  B,  Levitt,  M,  Linzer  and  B.F,  ^^ye^s 

SUMMARY 

Optical  reflectivity  measurements  of  shock  thickness  in  argon 
have  been  extended  to  much  higher  Mach  numbers.  Earlier  shock  thickness 
measurements  at  Mach  numbers  up  to  1.55  were  found  to  be  in  good 
agreement  with  Navier-Stokes  theory.  The  present  measurements  show 
that  this  agreement  becomes  progressively  less  good  as  the  Mach  number 
increases*  The  shock  thickness  in  argon  appears  to  reach  a  maximum  ab 
about  Mach  =  3  on  the  basis  of  preliminary  experiments  conducted  at 
M  =  4,99  where  the  obseived  thickness  agrees  with  the  bimodal  model 
of  Muckenfiiss. 

The  optical  reflectivity  method  has  also  been  applied  to  hydrogen- 
oxygen  detonations  and  strong  shock  fronts  in  pure  oxygen.  These 
measurements  suggest  that  rotational  relaxation,  which  is  known  to  be 
rapid  at  low  temperatures,  may  become  slow  at  high  temperatures.  It 
also  appears  that  in  detonation  waves  heating  effects  may  extend  into 
the  region  of  the  initial  compression  process. 


I .  Introduction 


The  optical  reflectivity  technique  has  made  it  possible  in  the  past 
to  study  the  thickness  of  shock  fronts  at  strengths  up  to 
M  =  1.55  and  in  lhat  region  it  was  found  that  in  argon  thicknesses  were  in 
accord  with  the  Kavler-Stokes  equations.  In  polyatomic  gases  it  was  pos¬ 
sible  to  measure  the  contribution  from  the  bulk  viscosity  and  hence  to 
determine  apparent  rotational  relaxation  times,  5*3  collisions  being  re¬ 
quired  to  atta'n  rotational  equilibration  in  the  case  of  nitrogen  and  fewer 
for  heavier  molecules.  In  the  present  studies  measurements  were  extended 
to  much  stronger  shocks  and  to  detonation  waves. 


II.  Experimental 

The  optical  reflectivity  method  is  based  on  the  fact  that  light  can  be 
reflected  from  a  shock  front  and  that  the  reflectivity  is  a  sensitive  func¬ 
tion  of  the  density  profile.  It  has  been  shown  (3)  that  when  R^'^1;  the 
reflectivity  of  a  shock  front  is  given  by  the  expression 

R(y)  =  I  F(y)  I 


where 


T 

,00 


F(y)  = 


-  ( 
j 

-'-00 


(1) 

(2) 


and  y  =  2  cosO/^  (5) 

Here  9  is  the  angle  of  incidence,  ^  the  wavelength  of  the  radiation  and  n 
the  index  of  refraction  as  a  function  of  the  distance  x  along  an  axis  perpen¬ 
dicular  to  the  shock  front.  The  index  of  refraction  is  related  to  the 
density  by  the  exi  ression 


/^(x)  _  n(x)-l 
yO  0  (no-1) 

where  the  o  subscript  refers  to  the  reference  state.  If  R(y)  could  be 


ih) 


measured  over  a  sufficient  range  of  y,  n(x)  and  hence^  (x)  could  be  de¬ 
termined;  in  practice  this  is  not  possible  and  an  analytic  form  with  adjustable 


2. 


constants  must  be  fitted  to  the  data. 

Ifcitil  now  we  have  used  the  function 

n(x)  =  n^  +  ^5) 

in  which  n^  is  the  refractive  index  of  the  unshocked  gas  and  L  is  a  parameter 
measuring  the  shock  thickness. 


L  = 


(6) 


For  this  model  the  reflectivity  is, 

R  = 


from  Eqs.  (l)  and  (2) 

2  (L/r^oo-  'VA)^ 

sill,?  V --'i.cosO/^  ) 


(7) 


The  function  log  R’  =  log  ,  \  plotted 

8inh^(7^^Lcos(?/>  ) 

in  Fig.  1  as  a  function  of  log  (X/I*  cos  c). 

This  model  (Eq.  5)  is  exact  frr  ?Ja*/ier-*Stokes  equations  with  constant 

coefficients  as  well  as  for  the  biuiodal  xhaory  of  Mott-Smith  and  Muckenfuss. 

It  is  a  limiting  case  of  the  expression  obtained  from  Sherman's  solution  of 
the  Burnett  equations,  (ll,12) 


^  5  +  ( 

5  -  ^  tanh(2x/L) 


where 


(8) 

(9) 


/'  15(M^-1) 

V  '!5+5M^ 

when  the  shock  strength  parameter  f-  is  small. 

The  technique  enployed  is  to  d  v’ ?rmin3  the  value  of  the  parameters  (/dn) 

and  L  in  Eq.  (7)  which  gives  the  best  fit  to  the  experimental  reflectivity 

function.  It  should  be  noted  that  although  L  is  a  "maximum  slope  thlckensp’’, 

the  method  does  not  itself  measure  this  quantity.  It  measures  the  best  fit 


of  a  function  of  the  form  of  Eq.  (5)  to  the  true  density  function. 

The  optical  system  employed  in  the  present  experiments  is  shown  in  Fig.  2. 


Radiation  from  an  Osram  HB0.109  super  pressure  mercury  arc  lamp  is  col¬ 
limated  and  iD5)inges  on  the  shock  front  as  it  passes  the  measuring  station. 


A  small  fraction  of  the  incident  radiation,  of  the  order  of  IC’®  to  10 
is  reflected  and  passes  successively  through  (a)  the  system  iris  which 
controls  the  area  of  the  shock  front  seen  by  the  detector,  (b)  a  focussing 
lens,  (o)  an  Interference  filter  with  a  bandwidth  of  about  7  and  a  peak 
wavelength  of  and  (d)  a  slit  which  controls  the  divergence  of  the 

beam.  It  is  finally  brought  to  focus  near  the  photomultiplier  tube. 

The  unreflected  portion  of  the  light  passes  through  a  similar  interference 
filter  and  is  received  by  a  phototube  whose  output  serves  to  monitor  the 


intensity  of  the  mercury  lamp.  In  most  of  the  experiments  both 


^  and 


were  kept  constant  and  the  reflectivity  curve  was  scanned  by  varying  the 
initial  pressure,  P^,  and  hence  L  which  is  related  to  it  by  the  scaling 


LiPi  =  LqPo  (■ 

The  apparatus  was  calibrated  by  measuring  the  reflection  from  weak 
shock  waves  in  argon  for  which  (/^n)  was  known  from  the  e:q)re8sion 


(I: 


and  L  was  known  both  from  theory  (5)  and  previous  measurements.  Huvever 
the  use  of  large  values  of  ^/L  cos  6  (Fig.  l)  made  the  calibration  re¬ 
flectivity  insensitive  to  the  value  of  L  used. 

The  accuracy  and  sensitivity  of  the  method  has  been  considerably 
improved  over  earlier  measurements  as  a  result  of  using  the  high  intensity 
source,  narrow  band  filters  and  of  refining  the  optical  system. 


III.  Shock  Front  Thicknesses  in  Argon 

The  results  of  recent  measurements  of  the  shock  front  thickness  in 
Ar(8)  sure  listed  in  Table  I.  Here  estimated  standard  deviation 


k. 


of  Lq^  thicituess  at  an  iaitial  pressure  of  one  atuiosphere,  1q  is  the  mean 
free  path  of  the  hsurd  sphere  gas  having  the  same  viscosity  as  Ar(lQ  =  0.69  x  10  ® 
cm  at  N.T.P.).  The  error  limits  are  the  90°/o  confidence  levels.  Fig.  3  shows 
a  ccanparlson  of  the  experimental  results  with  two  theoretical  curves.  In  the 
region  of  M  =  2.5,  the  shock  front  thickness  is  about  25°/o  greater  than  that 
predicted,  by  the  appropriate  Navier-Stokes  equation  but  agrees  with  the  bimodal 
calculation  of  Muckenfuss  (lO).  Also  at  M  =  4.99  where  data  from  preliminary 
measurements  are  reported,  the  experimental  shock  thickness  is  in  agreement  with 
that  predicted  by  the  bimodal  calculation. 

IV.  Detonation  Fronts 

According  to  the  accepted  model  of  detonations,  the  reaction  zone  is  pre¬ 
ceded  by  a  shock  wave  which  heats  and  compresses  tht.  reactant  gas.  Subsequently 
the  chanical  reaction  raises  the  temperature  and  consequently  reduces  the  density. 
For  the  study  of  the  initial  shock  wave  the  reflectivity  method  is  particularly 
suited  since  its  effective  time  resolution  is  of  the  order  of  lO'^^sec  whereas 
the  required  resolution  is  estimated  (4)  to  be  only  lo”®  sec.  (the  time  elapsing 
before  the  chemical  reaction  sets  in). 

Distinguishing  the  reflected  radiation  from  that  due  to  the  luminous  gases 
behind  the  suock  front  was  the  prime  difficulty  in  applyirg  the  optical  reflec¬ 
tivity  technique  to  the  study  of  detonation  fronts.  This  problem  was  overcome 
by  the  improvements  previously  mentioned,  as  well  as  by  measuring  radiation 
reflected  from  the  back  of  the  shock  front  so  that  the  viewing  system  did  not 
"see"  into  the  reaction  products.  Furthermore,  it  was  found  that  the  luminosity 
increased  as  the  first  power  of  the  initial  downstream  pressure  whereas  the  re¬ 
flectivity  increases  somewhat  faster  than  the  square.  Thus  the  interference  from 
gas  luminosity  cculd  be  reduced  by  working  at  higher  initial  downstream  pressures. 
In  practice  it  was  found  possible  to  make  measurements  at  initial  pressures  down 


5. 

to  10  p.a.i.a.  Finally,  the  luminosity  was  reduced  by  adding  diluent  gases 
which  reduced  the  temperature  behind  the  detonation  wave. 

The  first  detonation  experiments  were  conducted  with  a  HgOOs  gas  mixture 
(9).  The  in^iortant  results  of  this  work  were; 

1.  the  detonation  front  could  be  regarded  as  an  optical  surface  for 
which  the  angular  dispersion  was  less  than  O.05  reidlan, 

2.  the  detonation  front  thickness  (9xl0"^  cm  at  M=4,6)  differed 
very  slightly  from  shock  front  thicknesses  usually  found,  and 

3.  the  apparent  density  change  (  ^  =  2.6  )  was  only  ^k°/o  of  that 

ri 

expected  if  the  gas  achieved  equilibrium  behind  the  shock  front 
(  ^  =  3*5  )•  On  the  other  hand  it  was  not  far  from  that  ex- 
pected  if  there  was  no  rotational  contribution  to  the  heat 
capacity  of  the  oxygen  =  2.4  ),  i.e.,  if  rotational  re¬ 

laxation  were  slow. 

The  preceding  result  led  us  to  consider  the  following  possibilities: 

a.  that  there  was  little  rotational  excitation  of  O2  in  the 
shock  front,  and  that  about  20®/o  of  the  equilibrium  rotational 
energy  was  transferred  in  the  shock  front  for  which  the  tempera¬ 
ture  reached  was  about  1900°K.  Equilibration  would  therefore 
require  more  than  100  collisions  for  O2, 

b.  that  heat  flux  into  the  shock  front  reduced  the  density  change, 
perhaps  giving  a  profile  resembling  that  shown  in  Fig.  4  in  which 
case  the  theoretical  reflectivity  curve  of  Fig.  2  would  be  inap¬ 
plicable,  and 

c.  that  an  unlocated  experimental  error  was  responsible. 

To  check  these  possibilities,  experiments  were  conducted  with  the  gas  mix¬ 
ture  2H2:02JlOAr  for  which  rotational  relaxation  would  have  little  effect.  For 
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this  mixture, =  2. 50  with  no  rotational  heat  capacity  whereas  with  rota- 

A 

tional  equilibrium  for  Os  only,  ~  =  2.61.  The  results  obtained  are  shown  in 
Fig.  5«  It  is  clear  that  no  curve  drawn  through  these  points  can  be  superim¬ 
posed  on  the  theoretical  reflectivity  curve  of  Fig.  2.  These  results  are 
evidence  for  supporting  the  second  consideration  but  do  not  discredit  the  first. 
The  decrease  in  R’  at  the  highest  initial  pressures  suggests  that  the  density 
reaches  a  maximum  and  begins  to  fall  within  distances  comparable  to  the  thick¬ 
ness  of  the  shock  front. 

V.  Shock  Thickness  and  Rotational  Relaxation  in  Og 

As  a  further  check  on  the  hypothesis  that  rotational  relaxation  was  slow 
in  oxygen  at  high  temperatures,  preliminary  measurements  (9)  were  conducted  on 
shock  waves  at  M  =  3.1  in  oxygen  for  which  the  equilibrium  temperature  would  be 
834°K.  The  shocks  were  excited  by  detonating  acetylene-oxygen  mixtures  in  the 
driver  section  of  the  shock  tube.  These  measurements  yielded  an  apparent  shock 
front  thickness,  Lo,  of  approximately  3x10*®  cm  and  a  of  2.4+0. 2.  The 

latter  quantity  would  be  1.9  if  there  were  no  rotational  excitation  and  2.9  if 
rotational  equilibrium  were  achieved.  Consequently  there  appears  to  be  roughly 
70°/o  rotational  excitation  in  the  shock  front.  More  accurate  measurements  of 
rotational  relaxation  in  this  and  other  diatomic  gases  are  being  carried  out. 

VI.  Conclusions. 

It  has  now  been  possible  to  achieve  reasonable  precision  in  the  measure¬ 
ment  of  the  thickness  of  shock  fronts.  The  thickness  In  argon  reaches  a  mi ni mum 
at  about  M=3  and  at  M=4.99  Is  iu  agreement. with  that  predicted  by  the  bimodal 
theory  of  Muckenfuss.  Measurements  on  hydrogen-oxygen  detonations  as  well  as  on 
strong  shock  fronts  in  pure  oxygen  suggest  that  whereas  rotational  relaxation 
is  known  to  be  rapid  at  low  temperatures  it  may  become  slow  at  high  temperatures. 
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Finally,  it  appears  that  in  detonation  waves  heating  effects  may  e^rtend  into 
the  region  of  the  initial  compression  process. 
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Reciprocal  Shock  Thickness  (Mecm  Free 
Paths)  versus  Mach  Number 


FIGURE  3 


Mr 


14 


Experimental  Reflectivities 


